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We explore physics behind the horizon in eternal AdS Schwarzschild black holes. In 
dimension d > 3 , where the curvature grows large near the singularity, we find distinct but 
subtle signals of this singularity in the boundary CFT correlators. Building on previous 
work, we study correlation functions of operators on the two disjoint asymptotic boundaries 
of the spacetime by investigating the spacelike geodesies that join the boundaries. These 
dominate the correlators for large mass bulk fields. We show that the Penrose diagram 
for d > 3 is not square. As a result, the real geodesic connecting the two boundary points 
becomes almost null and bounces off the singularity at a finite boundary time t c ^ 0. 

If this geodesic were to dominate the correlator there would be a "light cone" singu- 
larity at t c . However, general properties of the boundary theory rule this out. In fact, we 
argue that the correlator is actually dominated by a complexified geodesic, whose proper- 
ties yield the large mass quasinormal mode frequencies previously found for this black hole. 
We find a branch cut in the correlator at small time (in the limit of large mass) , arising 
from coincidence of three geodesies. The t c singularity, a signal of the black hole singular- 
ity, occurs on a secondary sheet of the analytically continued correlator. Its properties are 
computationally accessible. 

The t c singularity persists to all orders in the 1/m expansion, for finite a', and to all 
orders in g s . Certain leading nonperturbative effects can also be studied. The behavior of 
these boundary theory quantities near t c gives, in principle, significant information about 
stringy and quantum behavior in the vicinity of the black hole singularity. 
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1. Introduction 



Some of the deepest mysteries in quantum gravity lie hidden behind horizons, includ- 
ing the nature of the spacelike singularities inside black holes. A better understanding 
of these may shed more light on cosmological singularities as well. Despite much recent 
work p],p|,|3|j^,|5|,P,[7|,|8|,|9|, p!0| , p!l|JT^ , p!3| , p!l| , these issues remain mysterious. More generally, the 



existence of horizons and resolving the puzzles of information loss and quantum decoher- 
ence require that the description of physics inside the horizon be linked somehow to the 
description outside. This connection-the central concept of black hole complementarity 
]T^,|16[-is still not well understood. 

As observations continue to bolster the case for cosmological inflation, this set of 
questions becomes more pressing. A complete description of eternal inflation |T7j probably 
requires an understanding of the many "universes" behind the de Sitter horizon of a single 



inflating patch |I||]||2(|. 



Despite the enormous progress in our understanding of quantum gravity in recent 
years, physics behind the horizon remains extremely puzzling. This is due in part to 



the holographic |T5| , |21[| nature of our most fully developed nonperturbative descriptions of 
quantum gravity, the brane bulk correspondences of Matrix Theory ]22[] and AdS/CFT 

These formulat 

ions describe physics from the point of view of an observer 



outside the horizon. For example, CFT correlators in AdS/CFT describe observables on 
the asymptotic spatial boundary of AdS space. In the thermal CFT that describes a large 
black hole in AdS/CFT f27j, these asymptotic observables apparently describe physics only 
outside the horizon. 

The brane bulk correspondences provide a fully self consistent quantum description 
of gravitational dynamics outside of the horizon, demonstrating that quantum coherence 
is not lost. But these descriptions appear to have little to say about behind the horizon 
phenomena. 

However, this is not the whole story. A number of authors have explored ways to 
extract behind the horizon information from boundary correlators P%| , |2!]| , |3U|pn| , |H2[] .B In 



particular Maldacena |3C§ argued that the boundary description of the eternal Anti-de 
Sitter Schwarzschild black hole consisting of one copy of the CFT on each of the two 
asymptotic boundaries could give some information about physics behind the horizon. 
In this approach, the one boundary thermal description is recovered by tracing over the 
Hilbert space of the other boundary CFT as in the real-time or thermofield formalism for 



thermal field theory [|4| . The thermal state counting entropy arises from the entanglement 



1 There is also large body of work addressing more general questions of extracting behind the 
horizon physics and holographic representation of horizons, as well as of generic local objects in 
the bulk; see e.g. H^HH,!!,^!!,!!®,!!!!!!!- 
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entropy of the pure entangled "Hartle-Hawking" state in which expectation values are 
computed. This realizes an old idea of Israel p5l . 



These ideas were explored in some detail for d = 3, the BTZ black hole, in [|3T]]. The 
finiteness of amplitudes (manifest from the boundary or bulk Euclidean point of view) was 
shown to arise in the bulk from ie regulation of the metric singularity inherent in analytic 
continuation from the Euclidean theory and from cancellations between future and past 
singularities. 

In addition, |3l]] gave a simple demonstration of how boundary correlation functions 
can probe physics behind the horizon by studying the correlator of two operators, one 
on each asymptotic boundary, each creating a large mass bulk particle. As the mass 
m — > oo, the correlator can be evaluated in the semiclassical geodesic approximation and 
is given by exp(— mC). Here £ is the proper length of the spacelike geodesic joining the 
boundary points (see Fig. 5a). Because the geodesic passes through spacetime regions inside 
the horizon, this boundary correlator reveals information about the geometry behind the 



horizon. This closely parallels a prior calculation done in the "geon" geometry |29 



As for all such quantities, an "outside the horizon" interpretation can be given. By 
tracing the Hartle-Hawking state, perturbed by the operator on one boundary, over the 
Hilbert space of the CFT on that boundary, a perturbed density matrix for the remaining 
CFT is created. This density matrix corresponds in the bulk to a modified boundary con- 
dition on the horizon, which can have nontrivial correlations with an asymptotic boundary 
field. The nature of this boundary condition is obscure from the "outside the horizon" 
point of view, but natural when physics behind the horizon is taken into account. 

The BTZ geodesies illustrated in Fig. 5a pass quite close to the singularity for moderate 
boundary times. But the correlation function is relatively structureless. Presumably this 
reflects the simple orbifold geometry of the BTZ black hole. The curvature of the BTZ 
geometry is constant except for a delta function at the singularity. 

Our goal in this paper is to explore further what information from behind the horizon 
can be found in boundary AdS/CFT correlators. To consider a situation with a more 
interesting geometry, we focus on AdS Schwarzschild black holes in d > 3 , and in particular 
in d = 5 where the boundary CFT is four dimensional Af = 4 super Yang Mills theory. The 
geometry near the singularity in such black holes approaches that of d > 3 Schwarzschild, 
and so the curvature diverges as one approaches the singularity. This affects the geodesies 
dramatically. 

We begin in Section 2 by exploring in detail the causal structure and spacelike 
geodesies in d = 5 AdS Schwarzschild. (Other dimensions greater than three behave 
similarly.) We find several unusual features. First, the Penrose diagram that encodes the 
causal structure of the geometry is not a square, unlike the d = 3 case. The lines describ- 
ing the singularity bend inward toward the center of the diagram. This allows a nearly 
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null geodesic beginning at the boundary at a finite nonzero time t c to "bounce" off the 
singularity and hit the other boundary at the symmetric point. This behavior, which does 
not occur for d = 3, is depicted in Fig. 5b. 

After submitting this paper we have learned that such features have previously been 
found in a general study of effectively two dimensional metrics [[16] @ . 

In Section 3 we turn to correlators in the geodesic approximation. If the bounce 
geodesic were to dominate, the correlator would become singular as t — > t c because the 
proper distance goes to zero as the geodesic becomes null. This would be a kind of "light 
cone singularity." But general considerations about the boundary field theory rule out 
this kind of behavior. In fact this geodesic does not dominate the correlator. There are 
in general multiple geodesies that connect the two boundary points (Fig. 12). At t = 
their proper distances coincide, creating a branch point in the correlator which behaves 
as £ 4//3 for small t. By studying various resolutions of this branch point, we show that as 
t increases from 0, the correlator defined by the boundary CFT is given by a symmetric 
sum of the two complex branches of this expression. Each of these can be attributed to a 
complex geodesic in complexified spacetime. At large t these complex branches reproduce 
the correct quasinormal modes for this black hole, which are complex, unlike the d = 3 
case. 

But the correlator is an analytic function of t and can be continued onto the real 
sheet. (This can be done in a computationally effective manner, as discussed in Appendix 
F.) On this sheet, the "light cone singularity" does appear. So the boundary correlator 
does contain information about the singularity, albeit in a subtle way. 

In Section 4 we extend the above analysis to include finite m, finite string length 
l s (a' ~ lg), and finite string coupling (g s ) effects. We take wrapped D-branes (giant 
gravitons P7| ) as explicit examples of heavy probe particles. The branch point at t = is 
smoothed out for any finite m, preventing analytic continuation to the bouncing geodesic. 
But at each order in the 1/m expansion the branch point persists and one can follow 
this geodesic and the accompanying fluctuation corrections into the bouncing domain. 
The 1/m corrections are given by a heat kernel expansion with coefficients related to the 
curvature and are large when the geodesic passes near the singularity. These quantities 
give a clear example of information about the singular geometry being coded into boundary 
CFT correlators. 

The main corrections at finite but small l s (with g s kept equal to zero) can be expressed 
as modifications of the supergravity field equations and hence of the background metric, 
dilaton, etc. fields. Such small modifications do not change the basic picture outlined 
above. The branch point at t = either persists or, more generically, splits into two 
nearby branch points. In either case we can follow the geodesic onto the real sheet. Of 

2 We are grateful to Jorma Louko for bringing this reference to our attention. 
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course, for t near t c , we expect the l s corrections to become large because the geodesic 
passes through regions of large curvature. 

At large l s the 't Hooft coupling A in the boundary gauge theory becomes small. In 
weak coupling perturbative gauge theory amplitudes are smooth at t = 0. There is no 
branch point. Therefore, for consistency, there must be a phase transition at finite A (for 
infinite m). 

Finite g s is a particularly difficult regime in which to study these phenomena. 
Wrapped D-brane masses are ~ l/<7s and so are finite here. This smooths out the branch 
point and would seem to prevent analytic continuation to the t c singularity. But by taking 
appropriate g s — > limits we can extract behavior around the singularity to all orders in 
g s , as well as study certain leading nonperturbative effects. 

These results demonstrate that a significant amount of information from behind the 
horizon, and in particular from near the singularity, is encoded in boundary theory corre- 
lators. In Section 5 we discuss the meaning of these results further. 

2. Classical geometry 

We will devote this section to the exploration of the classical Schwarzschild-AdS ge- 
ometry. As is well known, the AdS Schwarzschild spacetimes contain spacelike curvature 
singularities, event horizons, and timelike boundaries, all of which will play an important 
role in the discussion. First, we analyze the causal structure in detail, arriving at the sur- 
prising observation that, unlike the three dimensional case, for d > 3 the Penrose diagram 
of AdS Schwarzschild spacetimes cannot be drawn as a square, with both boundaries and 
singularities represented by straight lines. Next, we study spacelike geodesies in this ge- 
ometry, which, as we will discuss in the next section, provide a method for approximating 
the relevant boundary correlators. 

AdS Schwarzschild represents a d > 3 dimensional, two-parameter family of solutions 
characterized by the size of the black hole and the AdS radius. For simplicity we will 
focus on the large black hole (planar) limit, where the AdS radius is much smaller than 
the horizon radius. Also, we will from now on concentrate on d = 5. The geometry for 
finite-sized black holes, as well as in other dimensions d > 3, is qualitatively similar. We 
will contrast this with the quite different case of the d = 3 BTZ black hole ||31|| . 

The metric of the 5-dimensional AdS Schwarzschild spacetime is given by 



ds 2 = -f(r) dt 2 + 



dr 2 



+ r 2 dfl 2 



f(r) 



(2.1) 
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Fig. 1: Complexified coordinates for the AdS black hole: the time coordinate 
is complex in the extended spacetime, but with constant imaginary part in 
each wedge, as indicated above. The wedges are separated by the horizons 
at r = 1, where t diverges. Note that the real part of t increases upward in 
wedge I, downward in wedge III, and to the right (left) in wedge II (IV). 



where r + is the horizon radius and R is the AdS radius. We work in the limit of the 
infinitely massive black hole r + /R — > oo, where the metric simplifies. In particular, if we 
rescale the coordinates r — > -^r r, t — > £, measure lengths in AdS units so R = 1, and 
suppress the angular coordinates O that will not concern us here, the metric in the t-r 
plane becomes 

ds 2 = -f(r)dt 2 + </'•- 



f(r) = r 2 - 



M (2.2) 



There is a genuine curvature singularity at r = 0, and the boundary of the spacetime is 
approached as r — > oo. The Schwarzschild coordinates used in ( |2.2| ) also have a coordinate 
singularity at the horizon r = 1. This is not a big obstacle, since one can pass to Kruskal 
coordinates (T,X) which cover the full globally extended spacetime. (The Kruskal coordi- 
nate chart is constructed explicitly for this geometry in Appendix A). However, the global 
extension can be discussed more conveniently for our purposes by using four Schwarzschild 
coordinate patches (corresponding to the two asymptotic regions (I and III) with r > 1, 
plus the regions inside the black hole (II) and the white hole (IV) with < r < 1; cf. 
Fig.l). 

These four patches can be embedded in complexified Schwarzschild time, 

t = t L + it E (2.3) 
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where and tg denote times on the Lorentzian and Euclidean slices respectively. In each 
coordinate patch, t has a constant imaginary part, as indicated in Fig.l; we will define this 
to be for the right asymptotic region. Crossing a horizon shifts the imaginary part by 
i(3/A. Setting t = He produces Euclidean AdS Schwarzschild with a periodic imaginary 
time coordinate t^. The period of tE is f3, the inverse temperature T of the black hole. One 
can move a point from the boundary of region I to the symmetric point on the boundary of 
region III in Fig.l by rotating by half a period; that is, shifting t — > — t — i (3/2. The minus 
sign accounts for the opposite direction of time in region III. In the geometry (|2.2|) , /3 = it. 
The CFT correlators can be defined for complex time and are analytic, so such an extension 
is natural from the boundary side. The periodicity is just that of finite temperature field 
theory. 

2.1. Radial null geodesies 

We now turn to a major focus of this paper, the study of geodesies in the geometry 
(f2~2|). We consider a radial (O constant) path, x = (t(s),r(s)), where s is an affine 
parameter. Geodesic paths can be found by extremizing the action 



S = J ds x\ (2.4) 

where x = We can then find a unit speed parametrization by requiring x 2 = 0, +1, —1 
for null, spacelike, and timelike geodesies, respectively. 

Due to the symmetry described by the Killing field (-§i), there exists a conserved 
quantity E along each geodesic, 

E = tf(r). (2.5) 

To map out causal structure we first consider null geodesies. In a spherically symmetric 
spacetime, radial null rays will always be 45-degree lines in the Penrose diagram. Here 

■ 2 

x 2 = —f(r) t 2 + -jjpj = 0, which when combined with ( |2.5|) becomes 

r 2 = E 2 . (2.6) 

Let us consider the coordinate time t(r) along an ingoing radial null geodesic which 
starts from the boundary r = oo at t = to. From tfZlty and ( |2.6|) we have that 

<w = /^ / = < ° + f ( , 7) 

For geodesies which cross the horizon, integrating over the pole at r = 1 gives a constant 
imaginary part -if, so that for r < 1, 



2 



t(r) = to tan r tanh r — i — . (2. 
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(a) (b) (c) 

Fig. 2: Behaviour of radial null geodesies suggests the shape of the Penrose 
diagram (wavy lines represent the singularities, dashed lines the horizons, 
vertical lines/curves the boundaries, and the solid 45 degree lines the radial 
null geodesies): a) if geodesies meet at the singularity, the diagram can 
be drawn as a square; if they don't meet, b) with straight boundaries, the 
singularities are bowed in, or c) with straight singularities, the boundaries 
must be bowed out. 



The the last term is simply z/3/4, where again (3 = 1/T = n is the Bekenstein-Hawking 
temperature (cf. Fig.l) for this geometry. 

This result implies that a null geodesic which starts at the boundary at to = reaches 
the singularity r = at £(0) = f (1 — i). Since the real part of £(0) is positive, the geodesic 
hits the singularity off-center in the Penrose diagram; in fact it hits closer to the right 
boundary, where it originated. 

This immediately suggests the nature of the Penrose diagram, as illustrated in Fig. 2. 
If the Penrose diagram were a square, as in Fig. 2a, then radial null geodesies sent off 
from the boundary at £ = (defined as the horizontal line of symmetry of the Penrose 
diagram) would have to meet at the same point on the singularity. On the other hand, 
if the geodesies don't meet at or before reaching the singularity, then the singularity on 
the Penrose diagram must be drawn bowed in, assuming we draw the boundaries straight, 
as in Fig. 2b. Alternately, if the singularity were to be represented by a straight line, the 
boundaries would have to be bowed out, as in Fig. 2c, in order to reproduce the calculated 
behaviour of the geodesies. Note that the Penrose diagrams sketched in Fig. 2b and Fig. 2c 
are equivalent, since we can use a conformal transformation to straighten out the singular- 
ities of Fig. 2b, but at the expense of bowing out the boundaries. This will be discussed in 
greater detail in Appendix B, where the Penrose diagram for this spacetime is constructed 
explicitly. Such Penrose diagrams have previously been discussed in HH . 

Conversely, we can ask at what time does a null geodesic have to start from the 
boundary in order to reach the singularity at Re[£(r = 0)] = 0, i.e. in the "middle" of the 
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(a) (b) 

Fig. 3: Radial null geodesies on Kruskal diagram, a) for d = 3, and b) for 
d > 3 (radial null lines lie at 45 degrees, the horizons are represented by 
dashed lines, and the boundaries and singularities by the solid hyperbolas). 
In order for radial null geodesies to meet at the singularity at X = 0, they 
have to start at a) t = for the d = 3 case, and b) t = — t c < for the d > 3 
case. 



singularity. This will play an important role later. Denoting such a time by t c , we see that 

In general dimension d, one finds that t c = —tv (^(d — 1) tan (gzj/j j , which vanishes 
only for d = 3. Thus we see that, with the exception of the d = 3 BTZ black hole, all 
singularities in the Penrose diagrams of AdS Schwarzschild are bowed in as in Fig. 2b. 

In Kruskal coordinates (T, X), the future and past singularities of the black hole 
are the hyperbolas T 2 — X 2 = 1, while the AdS boundaries are given by T 2 — X 2 = 
—e 71 . (By contrast, in the d = 3 BTZ geometry, the singularities are at T 2 — X 2 = 1 
and the boundaries are at T 2 — X 2 = — 1.) The curved nature of the boundary of the 
Penrose diagram is reflected in these coordinates by the asymmetry in the radii of the 
two hyperbolas. This asymmetry allows a null geodesic starting at t c to "bounce off' the 
singularity. This is illustrated in Fig. 3, where the d = 3 case (Fig. 3a) is contrasted with 
the higher-dimensional case (Fig. 3b). In Appendices A and B we explain this in further 
detail, and construct the Penrose diagram explicitly. (The latter will be used in the next 
subsection to produce the Penrose diagrams in Fig. 5.) These issues are also discussed in 
4§ . 
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V(r) 




r 



Fig. 4: The "potential energy" in the mechanics problem for spacelike 
geodesies in AdS Schwarzschild in d = 3 and d = 5. The potential diverges 
at r = for all d > 3. 

2.2. Radial spacelike geodesies 

Let us now turn from the qualitative structure of the Penrose diagram to quantities 
which we can probe from the CFT side: in particular, the proper length along spacelike 
geodesies. In this subsection, we will consider purely radial spacelike geodesics.i The 

■ 2 

spacelike geodesies satisfy x 2 = ~f(r)t 2 + j^y = 1, giving 

E = tf(r), r 2 -f(r) = E 2 . (2.10) 

Here E, as before, is a conserved quantity resulting from the time translation symmetry 
of the geometry. We can think of it as the total energy of a particle moving in a potential 
V(r) = —f[r) (where the minus sign is due to the fact that we are considering spacelike 
geodesies). 

In Fig. 4, we plot the potential V(r) for d = 3 and d = 5 AdS Schwarzschild. In 
the BTZ case, the energy required to reach r = is finite. As a result, the proper time 
along the geodesic is large. By contrast, in the higher dimensional examples the potential 
diverges at the singularity, so that large E is required to approach r = 0. At large E the 
particle is moving fast and so very little proper time is covered. The result as E — ► oo 
is a geodesic that is null everywhere except in a vanishingly small region very near the 
singularity where it "bounces" off. As explained above, this occurs at to = t c . 

3 For completeness, we present in Appendix C an analogous treatment for spacelike geodesies 
carrying angular momentum. For the remainder of the paper, however, we need only the results 
about geodesies with zero angular momentum. 
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(a) Cb) 

Fig. 5: Symmetric spacelike radial geodesies of the AdS Schwarzschild black 
hole, plotted on the corresponding Penrose diagrams, a) d = 3: Because only 
finite E is required to come arbitrarily close to the singularity, none of these 
geodesies is nearly null, b) d > 3: the spacelike geodesies converge to a large 
E almost null geodesic which comes arbitrarily close to the singularity. 



The behaviour of radial spacelike symmetric geodesies (i.e. those which reach their 
minimum value of r at t = ^) in d = 3 and d = 5 is plotted on the corresponding Penrose 
diagrams in Fig. 5. As explained above, in the d = 3 BTZ case (Fig. 5a), the geodesies 
do not show any striking feature as the starting time to is varied. On the other hand, 
in higher dimensions such as d = 5 (Fig. 5b), spacelike symmetric geodesies exist only for 
a finite range of starting times, to £ ( — t c ,t c ). For starting times outside this range, no 
such geodesies exist. This is in accord with the expectation that the genuine curvature 
singularity of the higher dimensional black holes has a pronounced effect on geodesies in 
its vicinity, as opposed to the d = 3 BTZ orbifold singularity. 

Let us first consider the E = case. From ( |2.10| ), we see that t = 0; t(s) = is the 
only consistent solution. We have that r 2 = f(r). Such a geodesic crosses the Penrose 
diagram in the middle, so that it goes from r = oo at one boundary, down to the horizon 
r = 1, and back out to r = oo at the other boundary. Since it never penetrates inside the 
horizon, this geodesic can equally well be embedded in Euclidean AdS Schwarzschild. It 
is the geodesic that crosses over the tip of the cigar to the antipodal point = (3/2 on 
the Euclidean time circle (Fig. 8 ). 
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The proper length along such a geodesic is then given by 



max 



where r max — ► oo is the large r regulator standard in AdS/CFT |26| . Integrating this 
explicitly, we have 

£ = In (r^x + - l) « 2 lnr max + In 2 + 0(r^J (2.12) 

so that taking r max — > oo and subtracting off the universal logarithmically divergent piece, 
we obtain the regularized proper length 

£ reg = ln2 (2.13) 

In the future we will drop the subscript on C Teg . 

Let us now consider what happens as we increase the energy E. From (|2.10|) we see 
that the geodesic now penetrates a finite distance inside the horizon. Specifically, let r\ 
denote the smallest value of the radial coordinate reached along such a geodesic. This 
corresponds to the classical turning point of the particle motion, which is obtained by 
solving 

E 2 + f(r t ) =0 . (2.14) 

Since f(r) is positive outside the horizon and negative inside, < 1. For very high energy 
E ^> 1, the geodesic comes very close to the singularity, r{ <C 1, so that f(r) ~ — A, which 
means that rj ~ 4. 

In general, spacelike geodesies with arbitrary E and starting time to will not be sym- 
metric. We can specialize to this case, however, and obtain all others by t translation. 
For a given E, there exists a starting time to on the boundary in region I for which the 
geodesic is symmetric. In particular, for a symmetric geodesic the turning point r = ri 
has to occur at Re[t(rj)] = 0. The turning point must be in region II so Im[t(r^)] = —0/4. 
The geodesic ends in region III with coordinate t = —to — i/3/2. As mentioned before, the 
real part has a minus sign because the Schwarzschild coordinate time runs in the opposite 
directions in regions I and III. The geodesic equations give 

r°° t r°° e 

to-t(r t )=to + 0/A = - +dr = - - dr. (2.15) 

Jri r Jr % f(r) ^E 1 + f(r) 

For E = 0, we recover to = 0- One can also check that as E — > oo, to — > t c ; and in fact, 
for large energies one can see from examining the integral that 

(to-*c)~|. (2.16) 
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For each such geodesic there exists a mirror one obtained by taking t — > —t, E — > —E. 
The proper length, analogous to ( |2 . 1 1| ) , along a symmetric spacelike geodesic originat- 



ing at r max , going down to r i: and then back out to r max in the other asymptotic region, 
is 

£ = 2/ - (2.17) 

Upon evaluating and regularizing this integral, we obtain 

C = In [ . 2 ) . (2.18) 

As E goes to infinity the geodesic looks more and more null, and so at any fixed radial 
cutoff its proper length goes to zero. This is apparent already from Fig. 5b, where we can 
see explicitly that symmetric spacelike geodesies in the E — > oo limit approximate the 
null geodesic which bounces off the singularity. The regularized length along an E — > oo 
geodesic (which is roughly the difference between its proper length and that of the E = 
geodesic) goes to negative infinity at t c , 

£~21n(£-£ c ). (2.19) 
We will make use of this result in the next section, after we discuss the correlation functions. 



3. Correlation functions 

According to the AdS/CFT correspondence, AdS Schwarzschild in d = 5 is dual to 
M = 4 SYM theory at finite temperature. As before, we will consider the large black 
hole limit, so the CFT effectively lives in infinite volume. Standard CFT finite tempera- 
ture boundary correlators are dual to bulk correlation functions with insertions all in one 
asymptotic region (e.g., region I). There are other boundary correlators one might want to 
study. For example, one could put an operator on the boundary of region I, and another 
on the boundary of region III [pO| , |3T[] . We have discussed the geodesies connecting such 



points. In the field theory, such correlators have a natural representation in the thermofield 



double or real-time description of finite temperature field theory |44[ . 

In the real-time formalism, one takes the tensor product of two copies of the original 
field theory labeled by 1,2. The two copies are decoupled, and the total Hamiltonian is 

H tf = H®I-I®H*, (3.1) 

where H is the Hamiltonian for the original theory. We now construct the entangled state 

Z2 . 



12 



where \Ei, Ej) = \Ei) ® \Ej), and \Ei) are energy eigenstates. The state is a particular 
eigenvector of H t f with eigenvalue zero. Correlations between subsystems 1 and 2 are due 
to the entanglement in \if;). 

Operators which belong to subsystem 1 have the form A® I (where I is the identity 
operator), and will be denoted A\. Operators associated with subsystem 2 are defined in 
a similar manner, except with an additional rule of hermitian conjugation: 

A 2 = I®A^. (3.3) 



Standard thermal correlation functions may be written as an expectation value: 



(^\A 1 {0)B 1 {tM. (3.4) 

As can be easily seen from the form of l^), ( |3.4| ) is simply the thermal expectation value 
of A(0)B(t), evaluated in a thermal density matrix at inverse temperature f3. The state 
counting entropy observed in subsystem 1 is the entropy of entanglement of the state 
In field theory, no physical significance is usually attached to correlators involving both 
subsystems, but we can certainly define them; for example 



|;4i(0)fl a (tM. (3.5) 



In the finite temperature AdS/CFT correspondence, ( |3.5| ) has a simple interpretation [3C 
it corresponds in the bulk to a correlator between operators on the two disconnected 
boundaries of the spacetime. 

It is not hard to see that one can compute (|3.5|) by analytically continuing (|3.4|). 



MAi(0)fl 2 (*)M = MA^Bii-t - i/3/2)\4>)- (3-6) 

This is the analog of using complexified Schwarzschild time as discussed in Section 2. We 
see that these two sided correlators are just part of the information contained in ordinary 
thermal correlators, as a function of complex time. 

We can now use our knowledge of geodesies in AdS Schwarzschild to study correlation 
functions in the CFT. The prescription is to determine the 2-point CFT correlators via 
AdS/CFT from a computation of the bulk propagator. We assume at this point that we 
have a scalar field of mass m in the bulk, dual to some operator in the CFT whose 2-point 
function we want to compute. The bulk propagator is given by a sum over paths between 
two points in the bulk, with each path contributing e~ mC , C being the proper length of the 
path, suitably regulated. In the limit of large m, this sum will generically be dominated by 
the shortest geodesic connecting the two points. So the 2-point CFT correlator should go 
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like e~ mJ -' (up to ^ corrections), where C is the regularized length of the shortest geodesic 
between the two points^. 

In section 3.1 we show that a naive application of this prescription predicts a "light 
cone" pole in the 2-point opposite side correlator as t — > t c . However, simple arguments in 
the field theory show that this prediction cannot be correct, at least in the theory defined 
by analytic continuation from Euclidean space. In Section 3.2 we address this puzzle by 
examining the branch cut singularity in C at t = that is due to a coalescence of geodesies. 
We first resolve the branch cut by taking r max large but finite. We find that the unambigu- 
ously correct geodesic in Euclidean space bifurcates into two complex geodesies at small 
Lorentzian time t^. These geodesies must determine the correct answer. In section 3.3 we 
elucidate this further by moving this bifurcation onto the Euclidean section. This involves 
taking the black hole mass large but finite. Euclidean correlators are unambiguously deter- 
mined by the shortest geodesies at large m, and so we can determine the correct choice of 
geodesies through the bifurcation. The dominant geodesies continue to the complex ones 
mentioned above as one goes to the Lorentzian section. However, at large m it is possible 
to analytically continue correlators around these bifurcations by following geodesies that 
are no longer dominant. This allows us to find the t c singularity on a secondary sheet of 
the branched correlation function. This t c singularity gives information about the black 
hole singularity. It shows that such information is encoded in CFT correlators, albeit in a 
subtle way involving analytic continuation. In section 3.4 we introduce a simple model of 
the bifurcation given by an ordinary integral. The properties of this model can be deter- 
mined precisely and agree with the picture we have obtained. We demonstrate here that 
all orders fluctuation corrections around a given saddle point can be analytically continued 
as well. This indicates that additional information about the black hole singularity is also 
accessible. In section 3.6 we show that our results for the correlator agree with the known 
values for the quasinormal mode frequencies of the AdS Schwarzschild black hole. This 
gives independent confirmation of our analysis. 



3.1. A light cone singularity 

Recall from Equation ( |2.19|) that the regularized length of a symmetric geodesic with 
boundary time t diverges logarithmically at t c . Therefore the correlation function behaves 
like 

r 1 

e - mC rv, (3 7) 

We seem to have a pole of order 2m in the two-sided CFT correlator coming from the 
"almost null" geodesic connecting the boundary points. 



This approach and some of the pitfalls involved are discussed in 
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This immediately leads to a puzzle. We know that the real-time two-point correlation 
function, evaluated in a thermal state with inverse temperature /3, is 



(0(t)0H)>/3 = ^e-^"(E n |e- lH V(0)e lHt e lH V(0)e-^|i? n ) 



E p -f3E n -2it{E n -E m ) \i |2 ( 3-8 ) 
c |rnm| 5 



where <f) nm = (E n \(p(0)\E m ) and we have assumed (jy = <j). To get the two sided result, we 
can either use the real-time formalism and the state or continue one of the times from 
the real value t to t — i(3/2: 

<V# 2 (t)<M-W> = W + i/3/2)<K-t))fi = J2e- §{En+Em) - 2lt{En ~ Em) l^nml 2 - (3.9) 

n,m 

Note that the terms in the sum are real and positive for t = 0; therefore 

\(<h{t)M-t)h\ < (Mo)Mo))p (3.io) 

for any Hermitian operator (J). However, (<^ 2 (0)^1(0))^ can be computed unambiguously 
in Euclidean space, using the geodesic approximation. It is certainly finite (since the two 
points are on opposite sides of the Euclidean cigar). Thus we cannot have a singularity in 
the opposite side correlator at any time, including t = t c . 

3.2. Tracking geodesies 

We begin exploring this puzzle by computing the boundary time as a function of E 
from the integral ( |2.15| ). The computation is detailed in Appendix D; the result is 

■±E 2 + iE + l 





(3.11) 

( 9 \ 

£ = ln 

where formula ( |2 . 1 8| ) has been repeated for convenience. We have not restricted ourselves 
to symmetric geodesies here. The time t is the difference between tf and ti, the final and 
initial boundary Schwarzschild times; t = tf — ti + i/3/2. We have chosen the origin so 
that t = describes endpoints in regions I and III at zero time. For symmetric geodesies 
t = -2t . 

These equations determine C(t) and hence the large m correlator (<p<p)(t) = 
exp(— mC(t)). They are central to our analysis. 
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Fig. 6: Plot of £ as a functions of E, normalized by t(E — > oo), for various 
values of r max : the solid line has r max = 2, while the dashed line represents 

'"max 




Expanding ( |3.11| ) around E = we find 

£ ~ E z 



£ ~ -£ 4 (3.12) 



£ ~ — £3 

The branch point at £ = signals novel behavior. There are three branches for £ and 
hence for the correlator at each £. It will turn out the gauge theory chooses a symmetric 
combination of the two complex branches. The £ c singularity lies on the real branch. But 
knowing the correlator on the complex branch will allow us to study the real branch by 
analytic continuation. Information about the black hole singularity is encoded in the gauge 
theory in this manner. The analytic structure of the Riemann surface defined by (|3.11 ) is 
discussed in Appendix E. 
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To track these branches, it is useful to resolve the branch point. To do this we take 
r max finite. In the gauge theory this corresponds to a finite UV cutoff. Equations ( |3.12| ) 
become 

t ~ E 3 - aE 

(3 ' 13) 

Here a, b > are constants that go to zero as r max — > oo. At time t = there are now three 
real solutions for E and C (see Fig. 6). These correspond to the three geodesies shown in 
Fig.7. 

We can most precisely define gauge theory correlators by evaluating them in Euclidean 
time and then analytically continuing to general time. To compute the t = two sided 
correlator this way we smoothly shift one point of a coincident point Euclidean correlator 
by a half period, i(3/2. At t = we can do the geodesic computation entirely in the 
Euclidean AdS Schwarzschild geometry. The correct geodesic is clearly the E = solution 
of (|3.13[) which is described in ( |2.11| ). 

Now let us go to small Lorentzian time by increasing the real part of t by a small 
amount. The correct geodesic here is certainly the small deformation of the E = solution 
with E ~ t/a. This corresponds to following the central branch of the cubic curve in Fig. 6. 
But at a time t ~ a 3 / 2 this root of the cubic annihilates with another root at the local 
maximum shown in Fig. 6. These solutions then become complex. Complex values of E 
correspond to geodesies in the complexified AdS Schwarzschild spacetime, where both r 
and t are complex. These can be described as the solutions of the equations of motion of 
the mechanics problem with complex energy E. 

Therefore the gauge theory answer for the correlator, defined by continuation from 
Euclidean space, must become complex after a certain time. The complex branches of 
( |3. 11| ) do not contain the t c singularity, and so this analysis explains its absence in the 
gauge theory. On the complex branches the second derivative of the real part of C(t) is 
negative at t = 0. In other words, the correlation function computed by following this 
branch will decrease as we move away from t = 0, consistent with the general field theory 
result. 



3.3. Finite mass black hole 

To understand the pattern of geodesies better, we now describe how to move the 
bifurcation point to Euclidean time. The advantage of studying the Euclidean correlator 
is that, in the sum over paths representation, the weight factor is always positive and so 
we can reliably predict that the geodesies with the smallest £ dominate. To move the 
bifurcation we consider the finite mass AdS Schwarzschild geometry and compute t(E) for 
its geodesies. The result for small E takes the form (Appendix D): 



17 



Fig. 8: Geodesies of finite mass Euclidean AdS Schwarzschild with finite 
r max . The E = geodesic goes over the top and has larger C than the two 
geodesies with nonzero E. 



t~E 3 -aE 

(3 ' 14> 

Here a < 0, unlike the finite r max resolution described by ( |3.13| ). The parameters a and b 
go to zero as the black hole mass becomes infinite. The parameter b is still positive. 

If t — ti, + it E is Euclidean (pL = 0), then E is pure imaginary. Let E = —iE, with 
E real. Equation ( |3.14| ) becomes 

t E ~ E 3 - lol E 

(3.15) 

£~-E 4 -\b\E 2 . 1 ^ 

This cubic is also described by Fig. 6 with the horizontal axis now being E. At large E, 
the extension of ( |3.11[ ) gives t E ~ (3/2 describing almost coincident Euclidean points. As E 
is decreased, entering the domain of validity of ( |3.15| ), a single solution exists, describing 
a single geodesic. But at a particular value of t E ~ |a| 3 / 2 corresponding to the local 
maximum of the cubic in Fig. 6, a pair of geodesic solutions is created. From ( |3.15| ) and 
the sign of b it is clear that these new geodesies have larger L than the original one, so they 
are exponentially subdominant. This situation persists until t E = 0, when the geodesic 
endpoints are antipodal. Here the three solutions of ( |3.15| ) are E = 0, E = ±yjo[. The 
two geodesies with E = ±-v/|o| have the same £, which is less than that of the E = 
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geodesic. This situation is illustrated schematically in Fig. 8 where we have taken r max 
finite for clarity. It is clear from Fig. 8 that the E = geodesic is a local maximum of 
proper distance. 

At = the correct prescription to compute the correlator is to sum over the two 
dominant (E ^ 0) geodesies. We then can continue into the Lorentz section by varying th 
away from zero. To compute the correlator we follow both of the local minimum geodesies. 
From ( |3.14j) we see that E for this pair of geodesies will now have to have nonvanishing 
real and imaginary parts. If we label the two geodesies 1 and 2, we have E\ = E^ and 
L\ = The real parts of the proper distances are the same, so the magnitude of their 
contribution is the same. So our final formula for the two sided correlator as a function of 
Lorentzian time (or more generally, in a complexified neighborhood of Lorentzian time) is 
given by 

(</></>) (t) = e"™^ + e" m£2 . (3. 16) 

Because we have resolved the geodesic bifurcations we have been able to give an 
unambiguous calculation of this correlator. Taking the mass of the black hole to infinity 
we see that the correct prescription using Q3.11|) is to sum over both complex branches 
with equal weight. 

3.4- A bifurcation model 



As we have seen, at the branch point described in (|3.12| ), the geodesies coalesce. This 
means that stationary points of the path integral action become arbitrarily close in path 
space. This indicates that an eigenvalue of the fluctuation operator around such a geodesic 
must go to zero. There is a soft mode in the fluctuation spectrum around the branch point 
geodesic. More explicitly, the soft mode corresponds to infinitesimal motions of E away 
from zero; because t(E) ~ E 3 this means that the boundary points are fixed under such 
motions. In the space of paths, one can argue that the coordinate E is a good nonsingular 
coordinate for this mode near zero. For example, the central point of the path is deflected 
by a distance proportional to E. Since the regularized length behaves like C ~ i? 4 we see 
that the fluctuation eigenvalue goes like d 2 C/dE 2 ~ E 2 ~ t 2 / 3 . This zero mode at t = 
will create divergences in the 1/m expansion there. 

The fluctuation dynamics near t = should be dominated by this zero mode. So we 
should be able to get a reliable picture of the correlator in this region by truncating the 
path integral to this one fluctuating mode. Of course it will be non-gaussian, but from 
our knowledge of C(E) we should be able to determine its dynamics. Representing the 
nonsingular coordinate for the zero mode E by x = iE, and setting r = — it we can model 
the path integral for t ~ by the ordinary integral! 



This integral is a special case of Pearcey's Integral [48|. 
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Fig. 9: Integration contours for various values of t: left to right, r = 
1, e l7F / 8 , e l7r / 5 . At r = e 171 "/ 8 , the integral crosses a Stokes line and picks 
up another saddle. 



7(r) = / e" m ^ (x) rfx , (3.17) 

•/ — oo 

where the truncated action is given by 

V T = ^x 4 -rx. (3.18) 
To motivate this choice for V T , compute its stationary points. 



VS(x" - 



X s -T 



X« = T- 



(3.19) 



V(x a ) = -?r* 



This is the analog of ( |3.12|) .i For r real, (£ Euclidean), the integrand in ( |3.17| ) is real 



and positive on the real x integration contour. In such a situation the large m saddle point 
analysis is straightforward. We can then study deformations of t away from Euclidean. 

We now study this model for complex r. We can still use saddle point approximation 
for large m, but now we have to determine carefully which saddles lie on the integration 
contour. More precisely, for complex t we must deform the integration contour so that it 
becomes a steepest descent path for the real part of V T . Then we can evaluate I(r) by 
doing a saddle point expansion around all the saddles included on the path. A steepest 
descent path for the real part of V T is simply a contour of constant Im(V T ), so they are 
easy to identify. 



6 We will be able to study this integral in sufficient detail that we will not need to add the 
additional terms ~ x 2 to V that would resolve the cubic branch point. This can be done though. 
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Fig. 10: Positions of the Stokes and anti-Stokes (saddle exchange) lines in 
the complex r plane. The dashed lines are the Stokes lines, where the contour 
goes from crossing one saddle to two, and therefore picks up a sub-dominant 
contribution. The two saddles contribute equally on the wavy lines, so that 
the integral there has a kink (in the large m limit) as a function of r. 



For convenience, we restrict ourselves to Im(r) > 0, and write 

,3i6> 



r — re 



with 6 G [0, 4r], r > 0. The critical points are then 

x k =rz e + s 



(3.20) 



(3.21) 



where k = 0,1,2. Consider starting with real positive r, and then rotating r counter- 
clockwise into the complex plane. For 39 < ^, the steepest descent contour contains only 
the saddle point xq. For example, at 9 = 0, the steepest descent contour is simply the 
defining contour of the integral: the real axis. At 39 = |-, however, the imaginary parts 
of V T become equal at xq and x±, and the contour picks up the saddle point xi (Fig. 9). 
The line 3$ = n/8 — where the integral picks up a subdominant exponential in the saddle 
approximation — is known as a Stokes line |[49|| . This situation persists until 3$ = |7T, at 
which point the contour loses the point Xq, and only x\ remains (Fig. 10). 

By computing the real part of V T , we see that the saddle at xq dominates over x\ for 
3$ < 77, and vice versa for 39 > ^. The locus of this exchange of dominance is referred 
to as an anti-Stokes line. Right on the line both saddles contribute equally. This is the 
situation for the Lorentzian correlator discussed in the previous subsection. 

An interesting phenomenon occurs here. For | < 30 < f , the overall dominant 
saddle is £3, which is not on the integration contour, and does not contribute to the 
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integral in this approximation. This is similar to the phenomenon we encountered in 
the AdS Schwarzschild correlator, where the dominant branch (the t c geodesic) does not 
contribute. That we have a situation in which a dominant saddle does not contribute to 
the contour integral may at first seem surprising. In fact, it is easy to see that no matter 
how one tries to deform the contour so as to force it to go over this dominant saddle 
in a meaningful direction, one will pick up a larger extremum elsewhere on the contour. 
This extremum is of course only an extremum for V T restricted to the contour. It is not 
a critical point of V T as a function of complex r. However we do not have an analogous 
understanding of why the t c geodesic does not contribute to the Lorentzian functional 
integral. There must be other large contributions that cancel it. 

For m — > oo subdominant saddles make no contribution to this analytic function. 
The true m — > oo correlation function is a set of analytic regions, joined together in a 
non-smooth way at anti-Stokes lines when saddles exchange dominance. This is somewhat 
analogous to a first order transition in statistical mechanics. Computing the correlation 
function in one region and analytically continuing ignores these boundary lines, and is 
somewhat like continuing into a metastable phase. This analogy is imperfect, though. The 
bifurcation model illustrates that the t c saddle does not give a contribution like e~ mC to 
the path integral as would a metastable phase. 

Now consider I{r) for large but finite m. It is clear from the defining integral ( |3.17p 
that / is an entire function of r. The m = oo branch point at r = is smoothed out. 
The integral has an everywhere convergent power series and in a small disk |r| <C l/m*, 
I is roughly constant. When \t\ ~ 1 the saddle point approximation is very accurate for 
large m. Along the anti-Stokes (saddle exchange) line, the vertical axis in Fig. 10, there is 
a rapid but smooth transition from one saddle behavior to the other. The width of this 
transition is order 1/m. So there are two half plane regions where the large m saddle 
point behavior is very accurate, a thin joining region along the anti-Stokes line and a small 
roughly constant disk near the origin. These features smooth out the branch cut and so 
prevent us from following subdominant saddles by analytic continuation. 

A concrete method for performing analytic continuation is by computing power series 
coefficients around a point and then continuing the series variable. If we expand, for 
example, around a point away from the anti-Stokes (saddle exchange) line at large but 
finite m the contamination of the subdominant saddle point contribution that prevents 
continuation past the point of saddle exchange is of order e _m . Near the exchange line 
we can model the correction more accurately as e _m ( 1_r ) where the exchange line is at 
t = 1. The power series expansion of this correction goes like e~ m ^2(mr) k / kl . For 
fixed k the corrections are order e~ m but for k ~ m the corrections are order one. So the 

7 More precisely, we should consider the function L(t) = — — ln/(r) which has a smooth large 
m limit. 
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practical question about how large an m is required to be effectively infinite is equivalent 
to how many terms of a power series are required to accurately represent the infinite m 
function. The functions we are interested in have smooth large m limits so the number 
of terms required is m independent. This means that accurate analytic continuation can 
be done with results at large but finite m. Examples of this procedure for studying the t c 
singularity are discussed in Appendix F. 

There is an asymptotic series of 1/m corrections to 7(r) that can be computed by 
Feynman diagram techniques. Because of the soft mode, each of these terms is divergent 
at r = and has the structure l/(mr^) k . So each term lives on the same three sheeted 
Riemann surface as the leading answer. If we compute a given term around a dominant 
saddle and then analytically continue to a region where that saddle is no longer dominant 
we find the same answer we would get by perturbing around the subdominant saddle using 
a different integration contour. So the fluctuations we compute around the dominant saddle 
give us information about fluctuations around the subdominant one. We can clearly do 
the same thing with corrections in powers of a. We will use generalizations of this idea to 
argue that we can study stringy and quantum effects around the bouncing geodesic. 



3.5. The t c singularity 

We now have a reliable calculation of the boundary correlator at large m in the classical 
supergravity approximation as a function of complex boundary time t. The result is 
contained in Equation (|3.11|) , with additional information about which geodesies dominate. 



For Lorentzian time t = we sum over both complex solutions of (|3.11| ) with equal weight. 
The complex branches of ( |3.11 ) are free of additional singularities for all finite th- But 



the functions £i(t) are analytic in t. Given their values along any small segment of t 
it is possible to analytically continue them over the entire Riemann surface described by 
( |3.11| ) . The t c singularity can be found by such an analytic continuation from the dominant 
branch selected by the boundary CFT. In the discussion above of the bifurcation model 
we started from Euclidean t. Here let us give the argument starting from the Lorentzian t. 
Consider small t and hence use the expansion (|3.12| ), which gives us the relation E ~ . 
As we have argued previously, for small Lorentzian t the boundary CFT selects the branch 
E = e~ [* 3 1 . (We suppress the other complex branch which behaves similarly.) If we 
denote t = e %a t^ then we reach the branch with real t and real E (and hence real C) by 
taking a = 7r. On this branch there is a singularity at \t\ = t Cl due to the logarithmically 
diverging £ of the analytically continued geodesic, which "bounces" off the singularity. 
This t c singularity clearly reflects aspects of the black hole singularity. So we see that 
information about the black hole singularity is encoded in the boundary correlators, albeit 
in a subtle way involving analytic continuation. Despite its subtlety, analytic continuation 
can be computationally effective. In Appendix F we show how a modest number of terms 
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a) b) 

Fig. 11: Contours of constant imaginary part for the function t(E). There 
are three Im(t) = contours, which cross at the origin. The t c branch is the 
real axis. The other two give us the non-t c branches of the correlator. Pane 
b) is a magnified view of the singularity at 1 + showing only the Im(t) = 
contour. 



in a power series expansion of £ on a dominant branch is enough to extract precision 
information about the t c singularity. 



3.6. Quasinormal modes 



We can use the geodesic expression for the correlator derived above to compute the 
quasinormal modes of the AdS Schwarzschild black hole, in the large m limit. These agree 
with the known results ||50|1 , providing a check of our calculation, i 

To compute quasinormal modes we need the large t asymptotics of the correlator. We 
first note that along either complex branch E(t) must go to one of the fourth roots of —4 
as t goes to infinity. Specifically, let «i and a 2 be — 1 — i and 1 — % respectively; we can 
rewrite formula ( |3.11| ) as 



t 



I In 

4 



(E -a 1 )(E -o7i) 
(E - a 2 )(E - oT 2 )_ 



-iln 
4 



-(E - ai )(E - a 2 ) 
_ (E - at)(E - a 2 ) 



(3.22) 



We can expand around these various roots; for example, expanding E around a 2 gives 
t = c — (a 2 /4)ln(E — a 2 ), where c is a complex constant. Writing E — a 2 = re ld and 
requiring that t be real gives 

= hir + c', (3.23) 



We thank Chris Herzog for pointing out a mistake in a previous version of this section. 
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where d is a real constant. This demonstrates that the complex branches (where t goes 
to real infinity) are spirals in the complex E plane, converging on one of the roots (see 
Fig. 11). Coupled with the expansion C ~ — \h\(E — a) around any of the four roots, 
we have that as t goes to infinity, t ~ i(l + i)£ along one of the complex branches and 
t <~ 2 (1 — i)£ along the other. This gives the leading large m behavior of e~ mCl +e~ mC2 = 
e -m(i+i)t _|_ e -m(i-i)t This reproduces the first large m quasinormal modes of this black 
hole found in ||50|1 . Note that conventions differ between [|50| and this paper. 



To go beyond the first quasinormal mode, we argue as follows. Suppose we are on 
the complex branch that goes to c*2 as t goes to infinity. From ( [2.18| ) we then see that 
e~ c = (E — «2) 5 /(-E'), where / is analytic. We can invert and obtain 

E = ct2 + e~ 2C f(e~ 2C ) (3.24) 

with /(0) = 1 and / analytic. Substituting, we obtain t = i(l + i)C + g(e~ 2C ), with g 
analytic. Inverting once more we obtain 

C = (l-i)t + g(e- 2{1 - l)t ). (3.25) 

Now, the correlator is 

(<H>)(t) = e -Mm+d/m)f 1 (t)+-) ( 3 _ 26 ) 

where fi is the first 1/m correction. By extensivity and locality, we know that /i(t) = 
bot + h(t) where bo is a constant of order one determined by the fluctuation determinant 
around the geodesic (which we have not calculated) and h(t) is exponentially small in t 
because the world line path integral is massive in AdS Schwarzschild. Expanding gives a 
series 

oo 

(00) (t) ~ J2c n e- {rn+bo+2n ^ 1 - l)t +c.c. (3.27) 

This yields values for the quasinormal modes 

Wn = (m + b + 2n)(l±i) n = 0,1,2.... (3.28) 



These values agree with the large m quasinormal modes calculated in [SD[ if bo = — 1. 



The results in [EDI were obtained by a direct study of the radial equation for the AdS 



Schwarzschild geometry, the Heun differential equation. In some way the geodesic compu- 
tations summarized in ( |3.11j ) give a WKB solution to the Heun equation. 



4. Beyond the geodesic approximation 

In the previous section we studied the massive scalar correlator in the limit m — » oo 
where the geodesic approximation applies. We worked in the classical supergravity regime, 
a' = I 2 = and g s = 0. In this section we will discuss the more general situation. 
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4-1. Finite mass 

The scalar propagator has a path integral representation 



(<Kt 



_)0(f 2 )) = jf dT j Dx{s) exp ^-im/2 jf cte (i; 2 + l) j . (4.1) 



Here are points on the boundary and the integral is over paths in the bulk that connect 
these two points. For the radially separated points discussed in the previous sections 
(0(Ci)0(^2)) = {4 > 4 > ){^)- The mass m acts like 1/K, a saddle point parameter. Taking m 
finite means all paths are explored in the functional integral, not just the dominant saddles. 
So one of the key steps of the previous section-following a saddle point into a region where 
it is no longer dominant-will be problematic when m is finite. More explicitly, as discussed 
in section 6.5, the function L(m,t) = — — log(00)(t) defines a piecewise analytic function 
in the limit m —>■ 00 that can be continued past the anti- Stokes line where saddle point 
dominance is exchanged, but at finite m there is no clear way to follow the subdominant 
saddle past this line. Operationally, the "contamination" due to subdominant saddles must 
be controllably small to perform such analytic continuation. 

The coefficients of the 1/m expansion can be continued unambiguously and do contain 
significant physical information. From (|4.1f ) we see that this expansion is a small fluctuation 
or heat kernel expansion. There are two places where divergences are expected in the 
coefficients of (l/m) k . The first is t = 0, where multiple geodesies coincide. As discussed 
in section 3.3 we expect a soft mode in the fluctuations around the geodesic here. This 
will cause a divergence in the 1/m expansion. The terms in the 1/m expansion live on the 
same three sheeted Riemann surface as C. As discussed above, continuing these coefficients 
gives the small fluctuation expansion around the analytically continued saddle. 

The second place we expect divergences is on the bouncing geodesic for t ~ t c . Here 
the divergence is not due to a soft mode but to large nonlinearities present because the 
curvature is large near the singularity. We can expand the heat kernel ( |4.1| ) in a short time 
expansion 

(<Kvi)<Kv2)) ~Y^<* rW - > ( 4 - 2 ) 

fc=o \mj 

where r\i are bulk points near the singularity lying on the bouncing geodesic and Ry*> is 
shorthand for a curvature invariant that has scaling dimension (length) _2fc . Because the 
only length scale near the singularity is /, the shortest proper time between the geodesic 
and the singularity, we have ~ l~ 2k . The region of proper distance T where the 
geodesic is near the singularity is also ~ I. Putting in these estimates, (|4.2|) becomes 



k=0 ^ ' 
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This is reasonable since m I is the only dimensionless combination available near the sin- 
gularity. From ( |2.16| ) we see that r min ~ 1/E and E ~ l/(t — t c ). It follows immediately 
from the metric that / ~ r min ~ (t ~ tc) 2 - So we expect the 1/m expansion to have the 
structure 

<*>w~£*Uirn-F) ' (4 ' 4) 

We should be able to compute these coefficient functions on the primary sheet and the 
analytically continue them to the secondary sheet to study their singular behavior. This 
gives an example of recovering nontrivial information about the neighborhood of the sin- 
gularity from outside the horizon correlators. The information here is just the diverging 
curvature near the singularity, and it is clear from the above that the main contribution 
to the singular behavior comes from distances ~ / which can be made arbitrarily short (in 
supergravity approximation) by taking t — > t c . It is possible that the strengthening of the 
singularity with increasing k will have a recognizable signature in the correlator on the 
first sheet. 



4-2. Finite a' 

We now go beyond the supergravity approximation by letting the string length l s , 
( a' ~ I 2 . ) be finite. In the boundary SYM theory, this corresponds to taking the 't Hooft 
coupling A = gyM-N = (R/h) 4 finite. We still work in the classical limit ^ s — ^ 0, N — ► 00. 

When l s is finite, the behavior of correlators in the large m limit depends on which 
large m excitation we are studying. If we take a generic perturbative string state with 
many oscillator modes excited then the size of the excitation will diverge as m — > 00 and 
the correlator will cease to serve as a local probe. Instead we consider D-branes, which, at 



g s = are essentially pointlike, infinitely massive objects [51, 52,53] • In particular in the 
AdS$ x S 5 IIB string theory we study D3-branes wrapped around an S 3 submanifold of 
S 5 . Such a state is pointlike in the AdS§ space and realizes the BPS state carrying SO (6) 
charge, for large charge. These are the giant gravitons [l7j . 



The main effect of small but finite l s on such states is encapsulated in corrections to 
the supergravity Lagrangian. These include, for instance, corrections of the form (a') k R^ k \ 
The small parameter controlling the size of such terms is l s /l. For phenomena near the 
horizon I ~ R and this small parameter is of order l s /R ~ 1/A 1 / 4 . Such terms have a 
number of effects on the dynamics, including small shifts of the metric and dilaton fields 
[SIJ. 

We can estimate the effect of such shifts on the geodesic dynamics by following the 
discussion in the previous section about the finite mass black hole. Just as in that situation, 
the most general small, smooth shift in the geodesic equations can be encapsulated by 
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deforming ( |3.11|) and ( |2.18|) into 



t = E — aE 

Here a, b are small parameters whose signs are undetermined. As in the discussion follow- 
ing fl3.13|) and (|3.12|) this structure guarantees that the bifurcation, needed for analytic 



continuation, does not disappear at small but finite l s . At most the cube root singularity 
at t = splits into two square root singularities. 

This persistence of the singularity at small t creates a puzzle. Explicit field theory 
calculation of the two point correlator at small A (i.e., large l s ) shows no sign of such a 
singularity. So it seems there must be an m = oo phase transition at finite A. Perhaps we 
have missed some effect at large A that removes the singularity. On the other hand it does 
not look difficult to construct metric and dilaton modifications that remove the singularity 
for a finite size deformation. 

The nature of l s corrections for t near t c gives information about stringy behavior 
near the black hole singularity. As discussed above, we expect corrections to be controlled 
by the dimensionless ratio l s /l = A 1 / 4 /(t — t c ) 2 . So there should be large corrections for 
t~t c . 

Of course, this is a string tree level effect and can in principle be computed by string 
world sheet techniques. The D-brane trajectory is encoded in a boundary state. The 
analytic continuation necessary to study the real geodesic requires boundary states with 
complex D-brane positions. Such states have been studied recently in other contexts 
||55|| ||56|| . Perhaps one can circumvent the analytic continuation by directly formulating 
the conformal field theory in Lorentzian space, and then selecting the boundary state 
corresponding to a D-brane following the real geodesic. 

4-3. Finite g s 

We now relax the last remaining constraint and consider finite g s . A problem imme- 
diately arises because the masses of D-branes are ~ 1 / g s and so the mass m of the probe 
particle cannot be taken infinitely large. More precisely if the S 3 radius on which the 
D3-brane is wrapped is ~ it!, then its mass m is mR ~ (m s i?) 4 /(7 s ~ X/g s ~ N. If the S 3 
radius is string scale then mR ~ m s R/g s ~ AT 1 / 4 . If the SO(6) charge is Q then mR = Q.i 
The stringy exclusion principle |)7| and the giant graviton analysis reviewed above show 
that stability requires Q < N. As discussed above, at finite m the correlation function 



9 This is only approximately true because supersymmetry is broken in the black hole back 
ground, and the BPS conditions no longer strictly apply. But the semiclassical analysis in 
demonstrates that this relation continues to be approximately true for large Q. 



28 



contribution from the dominant saddle is "contaminated" by e~ m corrections from sub- 
dominant saddles that will eventually become dominant on crossing an anti-Stokes line. 
We can certainly study correlators as a power series in g s since ^> e~ m ~ e ~ x /9s as 
g s — > 0. Because the emission and absorption of virtual perturbative quanta cannot signifi- 
cantly alter the D-brane trajectory, we should be able to analytically continue the geodesic 
through the anti-Stokes line. This will allow us, in principle, to study the nature of string 
perturbation theory in the vicinity of the singularity by continuing each order across the 
anti-Stokes lines to t — > t c . Again, these quantities should in principle be computable by 
studying higher genus string diagrams with the appropriate boundary state. Of course 
there will be 1/m corrections to these order amplitudes which can be numerically larger 
than higher order g s corrections. But these are fluctuations about the main saddle, and 
will each analytically continue to the desired real geodesic. They can be distinguished from 
g s effects by their different parametric dependence on A, or equivalently Q. 

In the supergravity approximation at least, we might expect divergences in loop am- 
plitudes when interaction vertices approach the singularity. This issue was discussed for 
d = 3 in On general grounds, applicable here as well, such divergences should not 

occur. The correlators are manifestly finite in Euclidean SYM. They are analytic and so 
can be singular at most on a set of complex codimension one. So a generic divergence 
is not possible. From the bulk point of view the potential divergences in d = 3 are not 
present because the analytic continuation from Euclidean space induces a shift r — > r + ie 
which regulates the singularity. Interaction points are integrated across the singularity 
into another region to implement this. Remaining imaginary parts are cancelled between 
past and future singularities. Although d > 3 is more complicated we think it reasonable 
to expect a similar mechanism here. 

But the elimination of divergences need not eliminate large finite contributions. Since 
the geodesic near t c is very close to one singularity it is natural to expect large finite 
contributions in quantum corrections. We think this is what is being computed by the 
boundary theory near t c . 

We would like, at least in principle, to be able to obtain nonperturbative information 
about the singularity from the nonperturbatively well defined boundary CFT. 

Operationally, the size of nonperturbative phenomena whose continuation we want to 
study must be parametrically larger at the point we calculate than the e~ m contamination 
from the subdominant saddles that would become dominant at an anti-Stokes line and 
prevent analytic continuation of the "metastable" phase. For example, suppose we wish to 
study D-brane effects, for instance the pair production of small D-branes. Such effects are 
typically of size e _1//Ss . Using the largest stable giant graviton as a probe, e~ m ~ e~ x ^ 9s . 
Since e~ l ^ 9s is parametrically larger than e~ x ^ 9 % it should be possible to reliably compute 
such processes and then analytically continue them to the real sheet. Processes that may 
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be impossible to compute include NS brane production, expected to be of order e~ l ^ 9s . It 
may be possible, though, to isolate these effects by studying the boundary theory as an 
analytic function of g s and A. 

Processes where the nonperturbative effects are dominant should be much easier to 
study. As a first example, consider D-instanton effects that are the leading contribution to 
certain anomalous processes involving the bulk axion, related to the boundary 9 parameter. 
The effect of 9 on giant graviton correlators should behave like e~ 1 ^ 9a {fo{t) + 0(e~ rn )). 
Here we have normalized to the 9 = answer. The function fo(t) can in principle be 
computed on the primary sheet with parametrically small error. Then it can be continued 
to the secondary, real sheet. If D-instanton effects are enhanced near the singularity then 
there should be a signal in fo(t) near t c . 

As a second example, consider a charge Q wrapped D-brane "particle" of mass mi? = 
Q. Take Q ~ N. Now imagine a process where this particle fragments into two other such 
particles of charges Q\ and Q2, with Qi + Q 2 = Q and Qi ~ Q/2. The black hole is not 
a supersymmetric background so there is no BPS condition that prevents this. But this 
process is certainly nonperturbative. By considering the amount of D-brane that must 
be created and destroyed we have provisionally estimated that the rate for this process is 
~ exp(— Both initial and final particles follow geodesies of the kind we have been 
discussing, assuming that the final boundary operators at placed at the same point. These 
geodesies should continue smoothly onto the secondary sheet. Schematically, we expect the 
amplitude for this process to look like e~®(fo(t) + 0(e~®)) where we have normalized to 
the two particle charge Q correlator. Again, fo(t) is in principle calculable on the primary 
sheet with parametrically small error. It can then be analytically continued to t c . 

The kind of quantum information that can be obtained results from taking g s — > in 
a certain way. If one focuses on the behavior of amplitudes as t — > t c then perhaps the 
information that can be extracted about the singularity is some kind of double scaling limit 
that compensates for a shrinking g s by examining processes that grow large for t — > t c . 



5. Discussion 

The analysis of the preceding sections shows that a significant amount of information 
from behind the horizon, and in particular from the region near the singularity, is encoded 



10 In this entire discussion we have not considered the D-brane wrapped on a sphere in AdSs 
discussed [p8,59|. Such a state gets light near the black hole singularity and so it should have a 



strong effect on the analytically continued dynamics. It may be that the nonperturbative tunneling 



to this configuration discussed in [58.59] is the dominant mechanism for fragmentation near the 



black hole singularity. This deserves more study. 
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into boundary CFT correlators. An initial question we must consider is how this informa- 
tion can emerge from the hot, thermal, apparently featureless horizon. Our understanding 
of this issue is limited, but it is clear that analyticity plays a central role. It is analyticity 
that lets us reliably follow a geodesic into a region where its contribution is exponentially 
subdominant. One might think that our results are merely a consequence of the analytic 
nature of the classical AdS Schwarzschild geometry, which allows one to compute the com- 
plete metric from a small sample of it far outside the horizon. But clearly more than this 
is involved. We have argued in the previous section that a large amount of stringy and 
quantum information should be accessible by analytic continuation as well. It seems that 
it is the analyticity of the full quantum theory, or at least a limiting part of it, that is the 
central issue. 

Analyticity is not a central notion in classical general relativity. For instance, in the 
collapse of a matter shell to form a black hole the metric is not analytic, and only one 
asymptotic region exists. But outside the shell the metric is just the standard black hole 
metric. Correlators will in general not be analytic in such a spacetime. It seems that 
our methods may not apply to such cases. Even with analytic initial data, the same side 
correlators in a collapse scenario seem to differ from those of the eternal black hole by 
exponentially small terms, even though the geometry behind the horizon differs markedly. 

However, analyticity does play a central role in quantum field theory. Multipoint 
correlation functions are analytic functions of the point locations. So quantities computed 
in AdS/CFT will generically be analytic. We can envision setting up a collapse scenario 
in AdS/CFT by applying a large number of boundary operators at a fixed time far in the 
past. We then let this matter evolve. From the boundary point of view, we excite many 
gluons and then let them thermalize. We have created one of the microstates counted in 
the black hole entropy. We then can insert two very massive particle boundary operators 
as probes. Their correlator will be analytic in their time separation, and should go over 
to the thermal ensemble in the limit of large energy, which corresponds to large black 
hole mass. But even away from this limit the probe correlator is precisely analytic. In 
principle, we can test for periodicity in imaginary time, remnants of the t = branch point 
and of the t c singularity. We can look for exponentially small one sided contributions that 
become significant when continued to opposite sides. It may well be difficult to disentangle 
these small terms from the e~ m contamination that obstructs analytic continuation. The 
different parametric dependence and ability to study different microstates might make it 
possible to distinguish the two. This question may be easier to study in d = 3 where the 
real behind the horizon geodesic dominates for all t. If the effective geometry behind the 
horizon is modified then it should make a large effect on this two sided correlator. 

These questions should be related to the issue of recovering information that has 
fallen into the black hole. The lack of information loss in eternal black holes is signalled 
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Fig. 12: Symmetric spacelike geodesies of AdS Schwarzschild. The numbers 
indicate how many such geodesies pass through each point in that region. 



by Poincare recurrences ||30|| . which occur on time scales t r exponentially longer than t c , 
t r ~ exp(A^ 2 )t c . These recurrences may correspond to large fluctuations in the classical 
geometry [|30| . These occur at times ~ N 2 t c . It would be very interesting to connect 
these ideas to those discussed in this paper, although the enormous difference in time 
scales makes this challenging. 

As mentioned in Section 3, tracking the bouncing geodesic through the anti-Stokes line 
is a bit like following a metastable phase in a statistical mechanical system, although this 
analogy is imprecise. A related phenomenon has already been discussed in the AdS/CFT 
correspondence. Gross and Ooguri |)(| considered the expectation value of two facing 
Wilson loops of radius R and separation L. In the limit of infinite string tension A — > oo 
this quantity is determined by the minimal area string world sheet spanning the loops. As 
L is increased this surface jumps in a discontinuous, "first order" transition. At large L 
the disconnected surface dominates. But one can in principle continue through this "anti- 
Stokes line" and study the metastable connected surface until it pinches off into a kind 
of singularity. It would be interesting to consider various kinds of fluctuation corrections 
in this situation and see which ones can be continued through the anti-Stokes line to the 
pinch off point. 

We have argued that the real geodesic (which becomes almost null at t = i c ), con- 
tributes to boundary correlators in a subtle subleading way that can be exposed by analytic 
continuation to a secondary sheet. We might ask whether it ever provides the dominant 
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contribution to a physical quantity. We are not sure of the answer to this question. But, if 
we temporarily put aside the well founded concerns about measuring local bulk correlation 
functions in quantum gravity, and just compute them in classical supergravity approx- 
imation, we can find some answers. In Fig. 12 we have displayed the regions for which 
a symmetric bulk two point function is dominated by various geodesies. Inside regions 
3, 2, and the central 1, the most likely scenario is that the geodesic that dominates the 
correlation function lies on the real branch (real E) and so will evolve into the nearly null 
geodesic. For the dominant geodesic to get within a proper time / of the singularity, the 
correlator points have to be separated by a proper distance ~ /. For these quantities the 
large stringy and quantum effects discussed in the previous section should be important 
factors in computing the leading result, not a subtle subleading contribution. But of course 
the role of such bulk quantities is a deeply confusing one. 

We now ask what this new information can tell us about the singularity. Unfortunately, 
little is known about boundary CFT correlators at strong coupling. In fact, our results in 
the bulk make many new predictions about them, as has often been the case in AdS/CFT. 
Still, it is a useful exercise to ask what could be learned if complete correlator information 
were available for the CFT. 

In previous sections we concentrated on the two point correlator. But there is no 
reason not to study more general correlators (an example of this was discussed in Section 
4). In particular, consider the correlator with K operators on the asymptotic boundary of 
region I and K' on the boundary of region III, with a variety of SO (6) charge assignments. 
Each operator corresponds to a large mass particle (wrapped D-brane), and can be inserted 
at different points on the AdS sphere. Schwarzschild time runs in opposite directions in the 
two asymptotic regions, so the observer in region I can prepare an in-state of K particles, 
send the beam through the horizon, and then the other observer can observe an out- 
state of K' particles, assuming that we have analytically continued so the bounce geodesic 
dominates. If the boundary times of these operators are all near t c , the geodesies are 
almost null and hence almost on shell. A "meta-observer" could in principle compare the 
data of these two observers and compute a "meta S-matrix." This quantity would seem 
to contain significant information about the fate of organized matter as it approaches the 
singularity. 

But from the point of view of a global observer using Kruskal time, these correlators 
correspond to particle-antiparticle annihilation near the singularity, and the correlator is 
a kind of vacuum persistence amplitude. Its interpretation is less clear. 

A question one typically asks about singularities in consistent physical theories is what 
resolves them. Here the black hole singularity seems deeply connected to the boundary t c 
singularity in the analytically continued CFT amplitudes. Perhaps this t c singularity is 
smoothed out in the full CFT with finite A and g s . Of course the difficulties described in 
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the previous section in extracting finite g s information make it possible that this question 
can only be addressed in some double scaling limit. 

The definition of gauge theory correlators by analytic continuation from the Euclidean 
thermal theory must determine some boundary conditions for bulk fields at the black hole 
singularity. Some work has been done on this question in the d = 3 case [|6T 1, but d > 3 
is technically more difficult. We should be able to address this issue by tying geodesic 
techniques to WKB approximations of the wave equation. This may be connected to the 
possibility of defining different black hole systems by starting with different initial states 
than the Hartle-Hawking state. 

One important tool for analyzing the boundary CFT is to go to small 't Hooft coupling 
A (at N = oo) where weak coupling perturbation theory is valid. But, as mentioned in 
Section 4.2, at weak coupling there is no sign of the branch cut signalling the metastable 
phase on which we concentrate. As noted above it is not difficult to imagine Z s corrections 
to the metric and other supergravity fields that will remove the branch point and bouncing 
geodesic. It would be very useful to have an understanding of this from the gauge theory 
side. 
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Appendix A. Kruskal coordinates 

We have motivated the non-square nature of the Penrose diagram by studying radial 
null geodesies. Since the Schwarzschild coordinates do not cover the whole spacetime, we 
have used complexified coordinates in order to consider the globally extended spacetime. 
However, for the purposes of understanding the global structure, it is more desirable to 
pass to real, globally-defined coordinates. The most convenient ones to use are the Kruskal 
coordinates, which we now present. This will allow us to see the behavior of null geodesies, 
without having to use complexified coordinates and to deal with coordinate singularities 
at the horizons. 

Many of the issues discussed in this and the next Appendix have previously been 



addressed in a general context in [46 
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Suppressing the angular directions, the metric ( p72| ) is 

dr 2 



ds 2 = -f(r)dt 2 + j^ (A.I) 



with / = r 2 — 4j-. Now define the Tortoise coordinate 



o Z(r') 



+ C (A.2) 



with C an integration constant we will choose later, so as to make the Kruskal coordinates 
everywhere real. Letting u = t — r*, v = t + r#, we have (is 2 = — / dudv. This is singular 
at the horizon since the determinant vanishes; but we can perform a further coordinate 
transformation U = —e~ 2u , V = e 2v which will cast the metric into an extendible form. 
In these new coordinates, with f(r) and r*(r) being implicitly functions of U and V, we 
have 



ds 2 = -je- 4r *dUdV (A.3) 

Letting U = T — X and V = T + X, we can write this in a more familiar form ds 2 = 
g(T,X) (—dT 2 + dX 2 ). Note that by our choice of coordinates U(u) and V(v), the factor 
g(T, X) in front of (—dT 2 +dX 2 ) is nonsingular at the horizon: although / by itself vanishes 
at the horizon, the e -4r * factor cancels this, so that the conformal factor g remains finite 
across the horizon. Therefore, we can extend (U, V) from U < 0, V > 0, to U y V £ 
(—cxi, oo), subject to the constraint < r < oo. 

We easily derive the transformation laws0 between the old and the new coordinates: 

-e 4r *W =T 2 -X 2 (A.4) 
T 

2t = tanh _1 — (A.5) 
X 

Note that r* acquires an imaginary part as the integral passes through the pole at 1. It 
is correspondingly convenient to define C = ^f-, which makes r* real outside the horizon. 
More specifically, we have 

,4r, _ ( r ~ 1 ) „2 tan" 1 r 
r + 1 



Using (|A.4|) , we can now read off the curves of constant X 2 — T 2 labeling a constant 



r surface. In particular, it is easy to see that the singularity (r = 0) corresponds to 
T 2 —X 2 = 1, while the boundary (r — > oo) of AdS Schwarzschild is given by X 2 — T 2 = e 7r . 



li 



For general functions f(r) in (A.l) with simple zeros at r = r+, we would have T — X z 

T_ 

X 



. e /'(r+)r.(r) and 2 tanh- 1 § = f'(r+)t. 
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(a) (b) 
Fig. 13: Kruskal diagrams for AdS Schwarzschild, with a) d = 3 and b) 
d > 3. Radial null curves lie at 45 degrees; the horizons are represented by 
the diagonal dashed lines. The thick (spacelike) hyperbolas correspond to the 
singularity, the timelike ones to the boundary. The crucial difference between 
three and higher dimensions is that in the former case, the hyperbolas are the 
same distance from the origin, which translates into the geometric property 
that any radial null curve bouncing around the diagram ends up at the same 
point it started at. In higher dimensions this doesn't happen. 

Note that each of these equations is satisfied by two disjoint curves, related to each other 
by X — > —X, T — > —T. The full spacetime is then bounded by the four hyperbolas in 
Kruskal coordinates, 

-e w < T 2 -X 2 = UV < 1 (A.6) 

where the lower bound corresponds to the asymptotic (timelike) boundary, and the upper 
bound to the (spacelike) singularity. This is to be contrasted with the corresponding 
situation for the d = 3 BTZ spacetime, where a similar analysis would yield — 1 < T 2 — 
X 2 = UV < 1. 

Fig. 13 illustrates the contrast between the Kruskal diagram in d = 3 and d > 3 
dimensions. For the 3 dimensional BTZ spacetime, sketched in Fig. 13a, the singularities 
and boundaries intersect the Kruskal coordinate axes at the same distance. On the other 
hand, in all higher dimensions, the Kruskal diagram is as shown in Fig. 13b, where the 
singularities come closer to the center of the diagram than the boundaries. As we just 
computed, for the d = 5 large black hole, the ratio of these distances is e 71 " > 1. 
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This immediately implies that there cannot be a radial null geodesic starting on the 
boundary at t = (which implies T = 0) and hitting the singularity in the middle (X = 0). 
A more invariant way to say this is by considering a radial null geodesic (or a sequence 
of geodesies) , which "bounce aroundll! the diagram" , as shown by the thin solid diagonal 
lines in Fig. 13. It is easy to see that if the hyperbolas intersect at the same distance, as 
in Fig. 13a, then such a sequence of null curves ends up at exactly the same point where 
it started; otherwise, if the hypebolas intersect at different distances as in Fig. 13b, the 
geodesies end at a different point. Furthermore, if the boundaries are further than the 
singularities, these curves intersect inside the spacetime. As we will discuss in Appendix 
B, this simple geometrical property proves, in a more invariant way, that the Penrose 
diagram of AdS Schwarzschild cannot be drawn as a square for d > 3, but can for d = 3. 

So far, we have presented our results only in the large black hole limit. One may 
wonder whether for black holes comparable to (or much smaller than) the AdS radius, this 
effect does not go away. In fact, as we now explain, it gets larger. The smallest value that 
the magnitude of the lower bound in ( |A.6| ) can have occurs in the large mass limit, where 
this goes to X 2 — T 2 = e K . For smaller black holes this value increases, and in fact gets 

R 

arbitrarily large in the small mass limit: the boundary is then given by X 2 — T 2 = e r + 
with r + and R denoting the black hole and AdS radius, respectively. As we discuss next, 
on a Penrose diagram, such as sketched in Fig. 2b, the smaller black holes would have 
the singularities more bowed in (using the same conformal rescaling), or equivalently the 
boundaries more bowed out in Fig. 2c. But this is exactly what we would expect if we 
naively cut-off the asymptotically flat Schwarzschild Penrose diagram at larger and larger 
distances R (effectively corresponding to smaller and smaller black holes). 

Appendix B. Penrose diagram 

We have already discussed the qualitative features of the causal structure of AdS 
Schwarzschild black holes, but let us now confirm this by finding the explicit Penrose 
diagrams for these spacetimes. This will be achieved by finding appropriate transformation 
V — > v(V) and U — > u(U) of the Kruskal coordinates (V,U) (this guarantees that radial 
null geodesies will be 45 degree lines), such that the boundaries of the spacetime now lie 
at finite coordinate distance. 

12 Specifically, starting from any point on the right boundary, consider a future-directed ingoing 
null geodesic; from where this hits the future singularity, draw the past-directed outgoing null 
curve which eventually hits the left boundary; then the past-directed ingoing null curve which 
hits the past singularity; and finally the future-directed outgoing null curve which comes back to 
the right boundary again. 
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Fig. 14: Most general form of Penrose diagram with timelike boundaries 
and spacelike singularities being straight (thick lines on the diagram), since 
the horizons (dashed 45 degree lines) must join the opposite vertices. A null 
ray bouncing around (solid thin line) must come back to the same point. 

For example, letting V = tan | and U = tan ^ compactifies the spacetime into a 
region inside v, u G (— tt, 7r). Letting v = r + p and u = r — p, we can now check that 
the singularities UV = 1 are the straight lines r = ±f,p G (— §, §)■ However, the 
boundaries, given by (tan ^-j^ tan^-^) = — e 71 ', in the allowed region -0, -u G (— 7r,7r), can 
be explicitly checked to be "bowed out". On the other hand, if we try to straighten out 
the boundary, e.g. by taking V = e% tan | and U = e% tan ~, so that UV = — e n is given 
by p = ±f,T G (— f , f ), we find that the singularity must be "bowed in". An exact plot 
(generated by the latter transformation) was presented in Fig. 5b. 

One may ask whether one can make a more clever choice of u and v so as to straighten 
out both the singularity and the boundary. After all, the usual lore that "curves can be 
straightened out by suitable conformal transformations" should apply here as well. But in 
fact, we see that the present problem is sufficiently constrained so that this cannot happen. 

We can make a general geometric argument for why the singularities and the bound- 
aries cannot be drawn as straight lines, with the whole spacetime filling some compact 
region on the diagram. The first ingredient is to note that the horizons must be 45 de- 
gree lines, and that the singularities and boundaries must "meet" at these.0 (Note that 

13 By definition, horizons are the boundaries of the past of infinity (so the boundary must touch 
the horizon on the Penrose diagram), and any causal curve which enters the black hole must end 
at the singularity (so similarly, the singularity must touch the horizon in a Penrose diagram); 
and finally, the horizons themselves cannot form part of the boundary of the spacetime, as by 
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this is the reason the Penrose diagram could not be e.g. a rectangle, which might naively 
seem consistent with the observed properties of null geodesies.) The most general way the 
corresponding Penrose diagram with straight singularities and boundaries could look is 
sketched in Fig. 14. But one can easily show by a series of similar-triangle arguments that 
a null curve bouncing around this diagram (analogously to that sketched in Fig. 13) must 
come back to the same point. This is in contradiction to the properties of null geodesies, as 
we have demonstrated above, so the singularities and boundaries cannot both be straight, 
if the spacetime is to be mapped into a compact region. 

The last requirement is very important: we can easily construct maps which make 
both singularities and the boundaries straight, but at the cost of the "center" (where 
the two horizons intersect) being mapped to infinity. One such example is to let U = u 
and V = 1/v. The diagram of the spacetime would then not be compact, and would 
not qualify as a Penrose diagram. In retrospect, it is clear that one can't straighten 
out any two curves by conformal transformation, without letting them go to infinity. For 
example, consider two distinct curves which have the same endpoints. Whatever conformal 
transformation we apply, the two curves will still have the same endpoints. But clearly, 
in the flat geometry there is only a single straight line between any two points, so if 
the conformal transformation straightens out both curves without making them go off to 
infinity, it must be singular (i.e. non-invertible) . 

We should make one cautionary remark: we have shown that, keeping the boundaries 
straight, the singularities must be bowed in. The amount of bowing, however, does depend 
on the conformal transformation used. In other words, the crucial geometrical property 
we used was that radial null geodesies bouncing around the diagram end up at a different 
point from which they started. However, by a suitable conformal transformation, we can 
bring this point to appear arbitrarily close to the starting point, effectively making the 
Penrose diagram look almost, but not quite, like a square. 



Appendix C. Spacelike geodesies with angular momentum 

In this paper, we have been using only radial geodesies in AdS Schwarzschild. For 
completeness, let us now consider geodesies to ones carrying some "angular" momentum.0 
In this case, the behaviour of spacelike geodesies is altered significantly. In addition to 
timelike and radial components, the geodesies have an angular component, described by 

definition they lie entirely in the interior. 

14 We will use notation appropriate to finite-sized (spherically symmetric) AdS Schwarzschild 
black hole; in the large black hole limit, the angular momentum becomes the linear momentum 
along one of the translationally invariant directions. 
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the constant of motion L = <pr 2 . The radial equation ( |2.10|) is now modified to 

L 2 ' 



f' = E A + f(r) (CI) 
which describes a 1-dimensional motion of a particle of energy E 2 in the effective potential 

V eS (r) = -r 2 +L 2 + ±-^ (C.2) 

We see that rather than being repelled by the spacelike singularity, the spacelike geodesies 
which penetrate sufficiently close to the singularity are now attracted, and therefore termi- 
nate at the singularity. Namely, a spacelike geodesic with energy E and angular momentum 
L starting at r = oo can propagate to r = ri, where r{ is now a function of both E and L, 
given by the solution of the equation 

E 2 + f(n) (1-^=0 (C.3) 

Note that for E = 0, for any L, = 1 as before (assuming L < 1); however for E > the 
behaviour is more complicated, and depends on L. For L > 1, we see that > 1 since 
otherwise the LHS of ( |U.3| ) would be strictly positive. This means that such geodesies can 
never cross the horizon, and in particular can never connect the two boundaries of AdS 
Schwarzschild. Thus they are not relevant for us, and we will henceforth consider only 
geodesies with L < 1. For L < 1, the geodesies must propagate past the horizon (otherwise 
the LHS of (|C.3| ) would again remain positive), and there are now two possibilities: for 



high enough energy, the geodesic terminates at the singularity, whereas for energies below 
this critical energy, E < E c , it reaches only down to > L > 0, and reemerges in the 
other asymptotic region. The dividing energy E c is of course given by the maximum value 
of the effective potential in ( |C.2| ), E 2 = f{ri) {^r — . Note that this is consistent with 



our previous picture, since E c — >• 00 as — >• 0, i.e. as L — > + . In the small L limit, this 
behaves as E c ~ — ^. 

Confining attention to L < 1 and E < E c , let us consider the time to a t which a 
spacelike geodesic of given E and L has to start in order to be symmetric. By the same 
arguments as before, this is given by 

t (E, L) +i(3/4 = - = dr (C.4) 

Jn(E,L) /(r) ^E 2 + f(r) (l-£) 

As before, to vanishes for E = 0, and becomes more negative as E increases. Thus, to 
see the possible range of to for various allowed values of E and L, it suffices to consider 
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the part of the integral around for energies close to the saturating value E c . Since 
the effective potential peaks at E c , and ri(E — > E c ) approaches the r- value at which V e s 
attains its maximum, we can approximate the effective potential near its maximum by 
V e s = —f(r) (l — ~ E 2 — fi 2 (r — ri) 2 , where the coefficient fi 2 is some complicated 
function of L and the spacetime parameters. Substituting this back into ( |C.4| ) , we see that 
for energies E ~ E c , the integrand becomes %^ry • Since /(r) is well-behaved at r = ri 
for finite L, we see that the integral ( |(J.4j ) has logarithmic divergence at r^. Thus, unlike in 
the L = case where to was bounded by £ c (cf. ( p.!6|) ), for L > this bound disappears, 
and to can have arbitrary values. 

Let us now consider the behavior of the proper length along these symmetric spacelike 
geodesies carrying angular momentum. This is given by the generalization of ( 2.17 ), namely 

C = 2 (C.5) 
Jr * ^E 2 + f(r) (1 - 

where r max is the upper cutoff and rj is given by solution of (|C.3|) . Since we already 
analyzed the behavior as E increases, for simplicity we will now consider the E = case, 
to see how increasing L affects the proper length C. For E = 0, r{ = 1 independently 
of L, so we clearly see that increasing L must increase the proper length £ (since the 
denominator of the integrand in (|C.5| ) decreases). In fact, as L — >• 1, it is easy to see that 
the length diverges logarithmically at the near-horizon region. 



Appendix D. Derivation of the geodesic time t(E) 

In this appendix we derive the formula for the boundary time to as a function of E 
for the case of the symmetric spacelike geodesies. We first look at the infinitely massive 
black hole. The integral to be done is (|2.15|) , 

f°° E 

-to = / , dr (D.l) 

J n f(r) ^E 2 + f{r) 



where is the turning point and f(r) = r 2 — \. Letting r 2 = u, 7 = \ 1 + \E^ and 



u = v 7 — \E 2 , we obtain 



-t 







4- 

2-E 12 , we obtain 

\Eudu f°° \E(vi -\E 2 )~fdv 



[u 2 



-l) v /( W +iE 2 ) 2 -(l + ii? 4 ) Jl (hv-±E 2 ) 2 -l) 7 Vv 2 -l 

(D.2) 

Now let v = coshx. The integral becomes 



1 r°° (7 cosh x- \E 2 )dx E f°° (2u - ^-) dx 

2 E Jo 7 2 v 2 - 1 E 2 v + ±£ 4 - 1 = 4^ J v 2 - ^-v + \E* - 1 (D " >) 
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Note that the numerator is the derivative of the denominator with respect to v. This 
makes the decomposition into partial fractions easy, and the integral becomes 



j£ roo / 1 1 \ 

-t = — dx( + (D.4) 

4 7 Jo \v-a v — pj 

where a = - 2 — ^ — and (3 = ^— — are the two roots of the denominator of the previous 
equation. The integrals can now be done by substituting w = e x and writing cosh a; in 
terms of w. The answer is 



E 
2 7 



^a — V a \ 1 — /i a J flp — Up \l — flp 



(D.5) 



where \i x = x + y/x 2 — 1 and v x = x — V x 2 — 1. Using the actual expressions for a and (3 
in terms of E, we can simplify the answer to the following form: 

1 (±E 2 -E + l\ i (-±E 2 + iE+l . 
- to= in M— - -In ^^ = ^__ (D.6) 



Here we have subtracted the imaginary piece i/3/4 corresponding to crossing one horizon. 
One can check that this gives t c for large E and that it goes like E 3 for small E. This 
shows that there are 3 branches corresponding to the positive t direction starting at on 
the complex E plane. 

Recall that the result ( |3. 1 1|) quoted in Section 3.2 pertained to the time difference 
between two boundaries, so that for symmetric geodesies, t = — 2 to (defined to be real for 
real energies E). 

The case of finite mass black hole can also be done exactly, though it is computationally 
more involved. In particular, following the procedure outlined above for the metric ( |2 . 2| ) , 
we obtain 



ME) = , f „ {r + In ' (1 + ^ R2 ~* ER r + + H 



2{2r\ + R 2 ) I ~ E2 y R A + AR 2 r 2 ~ 4r 4 



•% + R 2 In 



(1 - E 2 ) R 2 + 2 i E R yjrl + R 2 + 2r\ 
1 + E 2 ) 2 R 4 + 4R 2 r\+Ar\ 



(D.7) 



Note that for t = and E = ^E, in the ^- -> limit, -t (E) reduces to ( jCTED with 



r+ ' r + 

barred t and E 



/ r 2 _|_^2 

Expanding around E ~ oo yields t c = —jf ^ + — which can be easily verified by 
considering radial null geodesic as in (|2.7| ). 
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Expanding around E m yields 



to(E) 



(2r 2 + + R 2 ) 2 



R 5 



E + 



R 5 (2r 2 + - R 2 ) 2 
2,{2r 2 + + R 2 Y 



E 3 + 0(E 5 ). 



(D.8) 



In passing we note the intriguing fact that the cubic term flips sign when r+ = -^=, which 
is exactly the same as the transition point when the specific heat becomes negative — that 
is, at the Hawking-Page transition |)!|]. Although the physical significance is not clear, it is 
rather suggestive that these small black holes (which are no longer approximately thermal 
states in the CFT) have a vastly richer geodesic structure. 

The important point, used in Section 3.3, is that t(E) now has the form 



with a < 0, so that the cubic degeneracy is resolved for Euclidean time. 

Appendix E. Analytic structure of E(t) 

We have given above a formula for t(E) in the infinite mass black hole limit. This 
formula defines a Riemann surface over the t-plane, and it is amusing to examine its 
analytic structure. Note that because of formula (|2.18| ) the analytic structure of E(t) is 
closely related to that of the large m correlation function (modulo extra branch cuts from 
the log in (|2.18|) ). To start, we find the locations of all the branch points of E(t). These 
must correspond to values of E for which dt/dE vanishes, and a quick examination of 
dt/dE shows that it vanishes only at E = 0, oo. 

Now t is periodic with period i/3 = in. For E(t), we have first of all the coincident 
points pole at iri/2 (with E ~ f _^ 2 ); our previous analysis shows that we have a branch 
cut at t = (with E ~ t^), and the t c pole at n/2. In fact this pattern of branch cuts 
and poles repeats itself, with branch cuts corresponding to points where E = and 
thus (looking at the residues of the logs in ( |3.11|) ) occurring at t = (p + iq)n/2 with p + q 
even, and 1/t poles corresponding to points where E = oo and occurring at t = (p + iq)n/2 
with p + q odd. Of course, these poles and branch cuts only occur on certain sheets of the 
Riemann surface. It is not too difficult to understand the sheet structure; the only hard 
part is seeing what happens at the branch points. There we have 3 sheets. It turns out that 
two of those (the non-t c branches) do not have any singularities or branch cuts. Moving 
far in Lorentzian time along these branches corresponds to E spiraling in to one of the 
fourth roots of —4, and thus getting far away from E = 0, oo, where the poles and branch 
points occur. Of course, the non-t c branches are the physical ones, so in some sense all the 
interesting analytic structure is very nonphysical. One possible interpretation for the string 



t(E) ~ E 3 -aE 



(D.9) 
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of t c poles is that they correspond to geodesies bouncing off both the singularities and the 
boundaries of AdS Schwarzschild multiple times. In order to bounce off a singularity the 
geodesic has to be spacelike whereas to bounce off the boundary of AdS Schwarzschild it 
has to be timelike, so this idea is speculative. It is worth noting, however, that the poles on 
both boundaries predicted above are precisely the locations where the geodesic bounces. 

Appendix F. Computing properties of the t c singularity 

In the text we argued that certain information about the black hole singularity is 
encoded in boundary correlators. In particular there is information in a singularity at 
t = t c on a secondary sheet of the two point correlator, which at m = oo has a third 
order branch point at t = 0. In this appendix we will show that this singularity can be 
studied in a computationally effective way, assuming that precise numerical information 
about the boundary correlator as a function of Euclidean time on the primary sheet is 
available. Of course such information is not available (for large A) from direct gauge 
theory computation at this time. But we feel it is useful to show by example how the 
process of analytic continuation can be made concrete given well defined gauge theory 
results. 

Specifically we assume that the a number of coefficients in the Taylor expansion of 
the gauge theory correlator (</>^)(£) are known. L{t) = ^ ln{4>4>) (t) has a smooth large m 
limit so we consider its expansion. For convenience we expand around euclidean antipodal 
separation, t = even though L(t) is branched there. The third order character of this 
branch point is known from our theoretical analysis or would be clear from numerical data. 

2 

So we write our expansion in terms of x = fs. For Euclidean time separations the large m 
evaluation of L(t) from a real geodesic is certainly correct. So we can expand ( j2.18|) and 
( |3. 11| ) to obtain model "numerical" gauge theory data. Stopping at order x 20 is sufficient 
to illustrate how things work. 

L(x) = 0.693147 - 1.36284a; 2 - 0.265333 a; 4 - 0.106447 x 6 - 0.0454717a; 8 

- 0.0199415 x 10 - 0.0090257 x 12 - 0.00421708 x 14 - 0.00202004 a; 16 (F.l) 

- 0.000984367a; 18 - 0.000485493 x 20 + . . . 

The analytic structure discussed in Section 3.2 and Appendix E implies that the t c 

2 /3 

singularity of L on the secondary t sheet will appear at x c = t c = 1.35128345. This 
singularity has the leading behavior L(x) ~ 21og(x — x c ). 

To find this singularity numerically we use standard techniques. First we differentiate 
( F.l|) so that the leading singular behavior is a simple pole, L'(x) ~ 2/(x — x c ). We then fit 
( |F.1| ) to a ratio of two polynomials P and Q of degree M and N, L'(x) ~ Pm(x)/Qn(x). 
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This is called a Pade approximant [E3| . The convergence of such approximants is nonuni- 
form, so we select one more or less arbitrarily. 

The [M, N] = [7,12] Pade has a pole at x = 1.35128349 (compared to x c = 
1.35128345) with a residue of 2.0028 (compared to 2) . These results are probably anoma- 
lously good. The [9, 10] approximant has a pole value of 1.35137. 

So in this limiting case a modest amount of "numerical" data is sufficient to study 
the singularity. But it is important to ascertain how stable this situation is. In particular, 
as discussed in Section 4.3, finite g s limits our ability to take m arbitrarily large. We 
know that at finite m the branch cut at t = is erased and we can no longer analytically 
continue to a secondary sheet. As discussed in Section 3.5, a signal of this erasure order 
one corrections to the order k power series coefficient when k ~ m. Here we see that k = 20 
is adequate to determine L, so we just need m » 20. As also noted in section 3.5 there 
will be an asymptotic series of 1/m corrections that are branched and can be analytically 
continued to the t c singularity. These contain interesting physics, as explained in Section 
4. However to study the leading behavior they must be handled in some way, perhaps by 
extrapolation in 1/m. But this may not be done with complete precision so we investigate 
our ability to still see the t c singularity with some small 1/m correction remaining. In the 
x variable this correction has a pole at x = and a higher order pole at x = x c . The 
correct procedure would be to expand around a point to 7^ and Pade approximate both 
x = and x = x c singularities. 

Rather than go through this entire procedure we model the phenomenon by adding 
the series expansion of a pole at x = |, m ( x li/ 2 ) t° the derivative of ( |F.1|) . This pole is 
much closer to the expansion point x = than x c is so its power series coefficients grow 
much faster than those in ( |F.1| ). The ratio of the coefficients at order n is approximately 
^(x c /(l/2)) n ~ ^(2.7) n At n = 19 for m = 100 the 1/m correction is about 10 7 times 
larger than the main result! Nonetheless the Pade approximant can efficiently separate out 
these different poles. The [7, 12] approximant (with m = 100) yields a pole at x = 1.35169 
with residue 2.0086, as well as the pole at x = 1/2. 

As explained in Section 4.1 we expect that the high curvature around the black hole 
singularity will create a leading 1/m correction of form l/{m{t — t c ) 2 ). The series expan- 
sion of this correction will eventually dominate over the terms in ( [F.l|) . The ratio of terms 
of order n will be roughly — n 2 . There exist techniques to disentangle such confluent sin- 



gularities in a Pade approximant |k| . But we will content ourselves with showing that 
this modest rate of growth allows us to ignore 1/m corrections for m ~ 10 4 . The [7, 12] ap- 
proximant with this addition gives a pole at x = 1.35061 with residue 1.990. This analysis 
illustrates how the apparently rather abstract notions of the m — > oo limit and subsequent 
analytic continuation to a singularity on a secondary sheet can be made computationally 
explicit. Of course the input data, Euclidean gauge theory correlators, do not yet exist! 
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